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Abstract 

We have measured the ratio R = T(Kl — > jj)/T(Kl — > 7r 7r 7r ) using the KLOE 
detector. From a sample of ~ 10 9 0-mesons produced at DA<£NE, the Frascati 
0— factory, we select ~ 1.6 10 8 if^-mesons tagged by observing Ks — > tt + it~ follow- 
ing the reaction e + e~ — > <j> — > K^Kg. From this sample we select 27,375 Kl — > 77 
events and obtain R = (2.79 ± 0.02 stoi ± 0.02 syst ) x 1CT 3 . Using the world average 
value for BR(K L -> 7r 7r°7r°), we obtain BR(K L -» 77) = (5.89±0.07±0.08) x 10" 4 
where the second error is due to the uncertainty on the 7r 7r°7r° branching fraction. 

PACS: 

keywords: 
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1 Introduction and experimental setup 

The decays K$ — > 77 and Kl — > 77 provide interesting tests [1] of chi- 
ral perturbation theory, ChPT. The dominant contribution to the K s — > 77 
decay is 0(p 4 ) and can therefore be computed with reasonable accuracy in 
ChPT. The 0(p 4 ) term vanishes for Kl — > 77 in the SU(3) limit. However 
large 0{p & ) contributions mediated by pseudoscalar mesons [2] are expected 
for Kl — > 77 with values depending on the amount of singlet-octet mixing 
[3]. A precise measurement of the Kl — > 77 decay rate is also of interest in 
connection with the Kl — > fi + fi~ decay. In fact the absorptive part of the 
decay rate, T(K L — > fi + fi~) a bs, is proportional to T(K L — > 77). This con- 
strains the dispersive part, T(Kl — > H + H~)dis and eventually the possibility 
of determining the V t d parameter of the CKM matrix [1]. Measurements of 
T(K S -> 77)/r(fT 5 -> 7r°7T ) and T(K L -> 77 )/r(K L -> 7r 7r 7r ) have been 
recently published by the NA48 Collaboration [4]. We describe a new mea- 
surement of T(Kl — > / ~f~y)/r(KL — > 7r°7T 7r ) obtained with i^-mesons from 

— > KsKl decays at DA$NE, the Frascati 0— factory. 

In DA$NE the electron and positron beams have energy E = m^j '(2 cos 9) 
where = 12.5 mrad is half of the beam crossing angle. 0-mesons are produced 
with a cross section of ~ 3 fib and a momentum of 12.5 MeV/c toward the 
center of the rings. 

The center of mass energy, W, the position of the beam crossing point (x, y, z) 
and the <p momentum are determined by measuring Bhabha scattering events. 
In a typical run of integrated luminosity / £dt ~ 100 nb^ 1 , lasting about 30 
minutes, we have 5W = 40 keV, Sp^ = 30 keV/ c, 5x = 30 /im, and 5y = 30 /im. 

The detector consists of a large cylindrical drift chamber, DC [6], surrounded 
by a lead-scintillating fiber sampling calorimeter, EMC [7], both immersed 
in a solenoidal magnetic field of 0.52 T with the axis parallel to the beams. 
The DC tracking volume extends from 28.5 to 190.5 cm in radius and is 340 
cm in length. For charged particles the transverse momentum resolution is 
Spt/pt — 0.4% and vertices are reconstructed with a spatial resolution of ~ 
3 mm. The calorimeter is divided into a barrel and two endcaps and covers 
98% of the solid angle. Photon energies and arrival times are measured with 
resolutions a E /E = 0.057 /y/E (GeV) and a t = 54 ps/y/E (GeV) © 50 ps 
respectively. The photon entry points are determined with an accuracy ai ~ 

1 cm/ ^Je (GeV) along the fibers, and ~ 1 cm in the transverse directions. A 
photon is defined as a calorimeter cluster not associated to a charged particle, 
by requiring that the distance along the fibers between the cluster centroid 
and the impact point of the nearest extrapolated track be greater than 3cr;. 
Two small calorimeters, QCAL [8], made with lead and scintillating tiles are 
wrapped around the low-beta quadrupoles to complete the hermeticity. 
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The trigger [9] uses information from both the calorimeter and the drift cham- 
ber. The EMC trigger requires two local energy deposits above threshold 
(E > 50 MeV in the barrel, E > 150 MeV in the endcaps). Recognition 
and rejection of cosmic-ray events is also performed at the trigger level, check- 
ing for the presence of two energy deposits above 30 MeV in the outermost 
calorimeter plane. The DC trigger is based on the multiplicity and topology 
of the hits in the drift cells. The trigger has a large time spread with respect 
to the beam crossing time. It is however synchronized with the machine radio 
frequency divided by four, T sync = 10.85 ns, with an accuracy of 50 ps. Dur- 
ing the period of data taking the bunch crossing period at DA<3>NE was T = 
5.43 ns. The T of the bunch crossing producing an event is determined offline 
during the event reconstruction. 



2 Data analysis 

The 0-meson decays into K$Kl ~ 34% of the time. The production of a Kl is 
tagged by the observation of a Ks — > 7r + 7r~ decay. Kl — > 77 and Kl — > 7r°7r 7r° 
decay vertices are reconstructed along the direction opposite to the K$ in the 
rest frame and required to be inside a given fiducial volume, FV . We call 
R = N^/N^o^o^o the ratio of interest. The numerators and denominators are 
found from: 

_ N obs — N bgd 

^trig ' ^tag ' ^FV ' ^sel 

where N t> s and Nb g d are the numbers of observed events and estimated back- 
ground, e t rig, ttag, £fv and e se i are respectively the trigger efficiency, the tag- 
ging efficiency, the acceptance in the fiducial volume and the selection effi- 
ciency for the two decays. The efficiencies e tag and e trig are equal at the few 
per mil level and cancel in the ratio R. Background and selection efficiencies 
must be separately determined. 

For this analysis the drift chamber is used to measure the K s — > it + it~ decay 
and to determine the direction of the Kl, the calorimeter is used to measure 
the photon energies and impact points and to reconstruct the Kl decay vertex 
by time of flight. 

The data sample was collected during 2001 and 2002 for an integrated lumi- 
nosity of ~ 362 pb -1 corresponding to the production of ~ 10 9 0. Details of 
the analysis can be found in reference [10]. Kl — > 77 events have a very clear 
signature, being the only source of ~ 250 MeV photon pairs that balance the 
momentum of the observed Ks- This allows the use of very loose selection 
criteria. On the other hand Kl — > 7r°7r 7r°, the dominant neutral decay, is 
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characterised by a large multiplicity of lower energy photons. The final error 
on R is dominated by the error on the number of Kl — > 77 events. 

Before full event reconstruction, the data are passed through a filter to reject 
machine background and cosmic ray events. As discussed later, this filter has 
a modest impact on the events of interest for this analysis. 

K s — > 7r + 7r~ decays are selected with the following requirements: 

- two tracks with opposite charge that form a vertex with cylindrical coordi- 
nates r v < 4 cm, \z v \ < 8 cm, and no other tracks connected to the vertex; 

- Ks momentum 100 MeV/c < p~K s = p n + + Pn- < 120 MeV in the rest 
system, and it + it~ invariant mass 490 MeV < M^+^- < 505 MeV. 

The Ks — > 7r + 7r - decay provides an unbiased tag for the Kl when it decays 
into neutral particles and a good measurement of the Kl momentum, px L = 
P<j> — Pk s , where p<f, is the central value of the e + e~ momentum determined 
with Bhabha scattering events. The angular resolution on the K L direction 
is determined from Kl — > /T + 7r~7r° events by measuring the angle between 
Pk l and the line joining the vertex and the 7r + 7r~ reconstructed vertex. The 
widths of the angular distributions are 5<p = 1.5°, 59 = 1.8°. 

The position of the K L vertex for K L — * 77 and K L — > 7r°7r 7r° decays is 
measured using the photon arrival times on the EMC. Each photon defines a 
time of flight triangle shown in Fig. 1. The three sides are the Kl decay length, 
Lk\ the distance from the decay vertex to the calorimeter cluster centroid, 
L 7 ; and the distance from the cluster to the vertex, L. The equations to 
determine the unknowns L K and L 7 are: 

L 2 + L 2 K -2LL K cos6 = L 2 
Lr/Pk + £ 7 = ct 7 

where t 1 is the photon arrival time on the EMC, (3kc is the Kl velocity and 9 
is the angle between L and Lk- Only one of the two solutions is kinematically 
correct. The value of L K is obtained from the energy weighted average, L K = 
Y^i Ei ■ Lki/ J2i Ei where i is the photon index. 

The accuracy of this method is checked by comparing in the Kl — > n + n^n 
decays the position of the K L decay vertex measured both by timing with 
the calorimeter and, with a much better precision, by tracking with the drift 
chamber. 

The resolution function, <j(Lk), is determined with Kl — > 7r 7r°7r° and Kl — > 
77 events by the distribution of the residuals L K i —L K where L K is the average 
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obtained with all the photons but the i . We measure for the Kl — > 7r 7r 7r° 
sample a^o^(L K ) = 2.06 - (0.16 • 10~ 2 L K ) + (0.19 • 10~ 4 L\) (cm) and for 
the K L — > 77 sample a^ 7 (L K ) = 1.73 + 0.0033 (cm). 

The Ki FV is defined in cylindrical coordinates as 30 cm < r t < 170 cm, \z\ < 
140 cm. The fraction of i^-mesons decaying in the FV is (31.5 ± 0.1)%. 




Fig. 1. The time of flight triangle 

The identification of the bunch crossing that originated the event is crucial to 
locate the vertex in space. An error by one bunch crossing period results in a 
displacement of the Kl vertex of about 33 cm and decreases the probability 
of correctly associating the photon clusters. The bunch crossing is determined 
by identifying one of the two pions of the Ks decay and by measuring its track 
length, momentum and time of flight. Thus an error of one (or more) crossing 
periods can occur if there is an incorrect track-to-cluster association or the 
track parameters are poorly measured. 

To minimise the number of events with an incorrect bunch-crossing assign- 
ment, we perform a consistency check of the time of flight of the pions along 
their trajectory l n measured with the DC, tec — ^n/PnlwC with the corre- 
sponding cluster time measured by the calorimeter, tEMC- Requiring \toc — 
tEMc\ < 2 ns for at least one pion, the probability of correctly identifying 
the bunch crossing is (99.4 ± 0.1)%. This additional cut retains 96% of the 
original Ks — > tt + tt~ event sample. The probability of identifying the correct 
bunch crossing was measured with a sample of Kl — > 7r + 7r~7r° decays where 
the position of the -k + tt~ vertex, r nn , is reconstructed by tracking in the DC 
and the position of the two-photon vertex, r 77 , by timing with the EMC. The 
difference r nn — r 77 is used to isolate the events in which the bunch crossing 
is incorrectly determined. 
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3 K L -> 7r 7r°7T selection 

The K L — > 7T 7r°7r decay has a large branching fraction, 21%, and thus has 
very small background. Given the large statistics we retain only 1 out of 10 de- 
cays. The selection of Kl — > 7r 7r°7r° events requires at least three calorimeter 
clusters with the following properties 

- energy larger than 20 MeV; 

- distance from any other cluster larger than 40 cm; 

- no association to a charged track; 

- Lk in the fiducial volume and \Lk% — Lk\ < ^<j{Lk). 

The main sources of inefficiencies are: 1) geometrical acceptance; 2) cluster 
energy threshold; 3) merging of clusters; 4) accidental association to a charged 
track; 5) Dalitz decay of one or more 7r°'s. The effect of these inefficiencies is 
to modify the relative population for events with 3, 4, 5, 6, 7 and > 8, clusters 
without significant loss of efficiency. Monte Carlo simulation shows that the 
selection efficiency is e se i = (99.80 ± 0.01)%. 

A comparison between data and Monte Carlo of the relative populations and of 
the distribution of the total energy, E = J2i Ei, shows that only events with 3 
and 4 clusters are contaminated by background. This is due to Kl — > 7r + 7r~7r° 
decays where one or two charged pions produce a cluster not associated to 
a track and neither track is associated to the K L vertex or to K L — > 7r°7r° 
decays, possibly in coincidence with machine background particles (e 1 * 1 or 7) 
that shower in the QCAL and generate soft neutral particles. 

To reduce this background, for the 3-cluster population we further require at 
least two clusters in the barrel with at least one of them with energy E > 50 
MeV and for the 4-cluster population at least one cluster in the barrel with 
energy E > 50 MeV. The probability to have a cluster with E > 50 MeV 
has been evaluated using the 6-cluster events. The probability of having a 
given number of clusters in the barrel depends only on geometry and has 
been evaluated by Monte Carlo simulation. We obtain e3,£>50MeV = (81.9 ± 
0.1)% and e 4) £ >5 oMeV — (98.3 ± 0.1)%. Additionally, an event with 3 clusters 
is accepted only if an additional cluster is found in QCAL within a time 
window of 10 ns with respect to the Kl decay time. The probability of such 
an occurrence is e qca i = (52 ±2)%. 

The K L — > 7T + 7r~7r° background is rejected by imposing a veto (track veto) on 
the events with charged tracks not associated to the Kg decay and with the 
first hit in the drift chamber at a distance of less than 30 cm from the position 
of the Kl vertex. The track veto also rejects about 60% of the Kl — > 7r°7r 7r° 
events with Dalitz decays. 
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Fig. 2. Kl — > 7T 7T 7r selection: distribution of the total energy for events with 3, 4, 
5, 6, 7 and > 8 clusters. Dots are data, shaded histogram is Monte Carlo simulation 
for Kl — > all channels. Data and Monte Carlo histograms are normalized to the 
same number of entries. 



Fig. 2 shows the distribution of the total energy for events with different num- 
bers of clusters together with the results of the Monte Carlo simulation. The 
relative fraction of events is shown in Table 1. The difference between data and 
Monte Carlo simulation for events with > 5 clusters is due to split clusters. 
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The contamination from accidental clusters originated by machine background 
is negligible. The residual background contamination in events with 3 and 4 
clusters is evaluated by Monte Carlo simulation and amounts to (18.6 ± 1.0)% 
and (7.0 ± 0.2)% respectively. 



Number of clusters 


Data 


Monte Carlo 


3 


0.37 ± 0.02 % 


0.35 ± 0.04 % 


4 


7.2 ±0.1 % 


7.3 ±0.1 % 


5 


31.5 ±0.1 % 


32.2 ± 0.2 % 


6 


57.4 ±0.1 % 


58.4 ± 0.2 % 


7 


3.3 ±0.1 % 


1.7 ±0.1 % 


> 8 


0.1 % 


0.03 % 



Table 1 

Fraction of events with at least three neutral clusters connected to the Kl decay 
vertex. 

A subsample of events has been processed and analysed without passing 
through the initial filter. The fractional loss due to the filter is found to be 
less than 10 -3 . The trigger efficiency for K s — > it + it~,K l — > 7r°7r 7r° events 
was measured in two different ways. A detailed description of the methods is 
given in reference [9]. The first method uses only the data and the informa- 
tion provided by the combined EMC + DC trigger. In the second method the 
Monte Carlo is used to evaluate the correlation between the EMC and the DC 
trigger showing that the correlation factor is very small. The results obtained 
with the two methods, e trigl = (99.88 ± 0.04)%, e trig2 = (99.90 ± 0.03)%, are 
in good agreement. Since the two methods are independent, and the results 
consistent, we combined the two results. 

The number of events is: 

M N^ + Ni + N 5 + N 6 + N 7 + A> 8 + N Dalltz 

^downscale ' ^trig ' ^sel 

where A3 and N' 4 are corrected for the background subtraction and the addi- 
tional cuts quoted before. N£, a i itz , a small addition of 0.46% of the total count, 
is obtained using the Monte Carlo result that 21.5% of the Dalitz decays are 
included in A3 + A 4 events while 60% of them are rejected by the track veto. 
We find A^o^o = 9, 802, 200 x (1 ± 0.0010 stat ± 0.0016^). 

To check the uniformity of the Kl — > 7r 7r°7r° vertex reconstruction efficiency 
throughout the FV we have studied the proper time distribution. From a fit 
to the distribution we find tk l = 51.6 ns with a statistical error of 0.4 ns [10], 
in good agreement with the value reported in PDG [5], t Kl = (51.7 ±0.4) ns. 
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4 K L — > 77 selection 

Kl ^ 11 events are preselected by requiring at least two calorimeter clusters 
with energy E y > 100 MeV not associated to tracks. For the two most energetic 
clusters we require: 

- total energy E 12 = E yl + E l2 > 350 MeV; 

- angle between the photon momenta projected onto the plane normal to the 
K L direction, ip > 150°; 

- time difference smaller than 15 ns; 

- K L decay vertex in the fiducial volume and AL K = \L K i—L K2 \ < 4<7 7T (L^). 

The geometrical acceptance and the selection efficiency are evaluated by Monte 
Carlo simulation. The values of the efficiency are shown in table 2. 



preselection 


efficiency 


£ 7 > 100 MeV) 
E 12 > 350 MeV 
Y> > 150° 
AL K < 4a 
At < 15 ns 


(92.6 ± 0.1 stat ) % 
(99.88 ± 0.02 stat ) % 
(98.4 ± 0.1 stat ± 0A syst ) % 
(98.5 ± 0.1 stat ± 0.2 syst ) % 
(99.89 ± 0.02 stat ) % 


total 


(89.5 ± 0.2 stat ± 0A syst ) % 



Table 2 

Efficiency of the Kl — > 77 preselection. 

With these cuts we obtain 1.7 x 10 5 events with a large background due to 
K L — > 7r 7r°7r° and — > 7r°7r° decays, K L — > 777 being negligible [11]. 

The signal is further selected using the two body Kl — > 77 decay kinematics. 
In fact, photon energies can be computed with better accuracy from cluster 
and decay vertex coordinates. The laboratory energy is obtained by boosting 
from the center of mass where E 1 = M K /2 to the laboratory. If p ri are unit 
vectors from the K L decay vertex to the cluster centroids, the photon energies 
are 

E > M K /2 
11 1k(1 - Pk P~,i ■ Pk) 

where (3k and px are computed from px L = P<f> — Pk s ■ The has energy 
E' = E' lX + E'^ 2 , and momentum p 77 = E' yl ■ p yl + E' l2 ■ p l2 . Fig. 3 shows 
the distribution of E' and of the angle a between p 77 and px L , together with 
the results of the Monte Carlo simulation. The data are fitted with a linear 
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combination of the Monte Carlo distributions for signal and background. The 
fit gives the relative normalisation for the two populations. 




E (MeV) a (degrees) 



Fig. 3. Kl — ► 77 selection: distributions of laboratory total energy E' (left) and the 
angle a between p 77 and px L (right). Dots are data, shaded histogram is Monte 
Carlo simulation for the signal, dashed histogram is Monte Carlo simulation for 
background and solid line histogram is the Monte Carlo simulation for signal and 
background. 

In order to reduce background we further require: 

- \E' — fj,'\ < 5a' where // = 510.0 MeV and o' = 1.8 MeV are evaluated from 
a fit to the E' distribution; 

- a < 15°. 

To extract the signal we fit the invariant mass M 77 distribution obtained using 
calorimeter cluster energies with a linear combination of the Monte Carlo 
distributions for signal and background. The result of the fit gives the number 
of events for the two populations. Fig. 4 shows the distribution of M 77 before 
and after the E' and a cuts. From a fit to the second distribution we find 
22, 185 ± 170 K L -> 77 events. 

The efficiency of the selection cuts are evaluated from the data using a sample 
of Kl — > 77 events with high purity (S/B ~ 10 3 ) selected by applying hard, 
uncorrelated cuts on other kinematic variables [10]. 

The systematic error associated with the selection cuts on (£", a) is evaluated 
by moving the cuts around the chosen values and fitting the invariant mass 
distribution. The maximum displacement of the measured value for the num- 
ber of signal counts is ±0.2% and ±0.3% for the E' and the a distribution 
respectively. The systematic error due to the background contamination has 
been evaluated by changing the shape of the background distribution used as 
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Fig. 4. Kl — > 77 selection: distributions of the invariant mass, M 77 , before (left) 
and after (right) the E' and a cuts. Dots are data, shaded histogram is Monte 
Carlo simulation for the signal, dashed histogram is Monte Carlo simulation for 
background and solid line is the Monte Carlo simulation for signal and background. 

input of the fit. The effect on the signal is 10 times smaller and produces a 
systematic error of ±0.3%. 



The filter and the trigger efficiencies are evaluated as for the analysis of Kl — > 
7r 7r°7T decay. The results are e filter = (99.93 ±0.01)%, e trig = (99.44 ±0.04)% 
where the statistical and systematic errors are combined in quadrature. The ef- 
ficiencies associated with the various analysis steps are summarised in Table 3. 



selection 


efficiency 


trigger 
filter 
preselection 
\E' -fi'\ < 5a' 
a < 15° 


(99.44 ± 0.04) % 
(99.93 ± 0.01) % 
(89.5 ± 0.2 stat ± 0A syst ) % 
(98.5 ± 0.2 stat ± 0.2 syst ) % 
(92.5 ± 0.3 stat ± 0.3 syst ) % 


total 


(81.0 ± 0.3 stat ± 0.5^) % 



Table 3 

Efficiencies for the selection of Kl — > 77 events. 



The number of events is A^ 77 = 27, 375 x (1 ± 0.007Q stat ± 0.0081 SJ/St ). For the 
ratio we find: 



R 



r(K L -> 77) 

T(K L -> vr 7r 7r°) 



(2.793 ± 0.022 stai ± 0.024^) x 10" 
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in good agreement with the recent result from the NA48 Collaboration T(Kl — > 
H)/V(K L -> 7r°7T 7r ) = (2.81 ± 0.01 stai ± 0.02 SJ/st ) x 1(T 3 [4]. 

Using the known value for the K L — > 7r 7r 7r° branching fraction, we obtain 
BR(K L -> 77) = (5.89±0.07±0.08) x 10 -4 where the first error represents the 
statistical and systematic error on R combined in quadrature and the second 
is due to the uncertainty in the branching fraction. A decay width of 

T(Kl — > 77) = (7.5 ±0.1) x 10~ 12 eV is in agreement with 0{p & ) predictions 
of ChPT provided the value of the pseudoscalar mixing angle is close to our 
recent measurement of 9 P = (-12. 9^1)° [12]. 
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